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Abstract 

In  this  article  we  discuss  successive  approximation 
method  for  a  system  of  random  Volterra  integral  equations. 

An  example  is  presented  to  implement  the  theory.  Kolmogorov- 
Smirnov  test  is  used  to  fit  a  distribution  of  the  solutions. 

1.  Introduction 

The  study  of  random  Volterra  integral  equations  and 
their  applications  play  an  important  role  in  the  area  of 
probabilistic  analysis.  One  of  the  main  reasons  is  that 
integral  equations  are  suitable  for  numerical  treatment. 

For  a  recent  survey  of  approximate  solution  of  random  inte¬ 
gral  equations  we  refer  to  Bharucha-Reid  and  Christensen 
[41  .  For  the  numerical  treatment  of  random  integral  equa¬ 
tions  we  refer  to  Bharucha-Reid  [3],  Becus  [2],  Christensen 
and  Bharucha-Reid  [5],  Lax  [8-10],  Tsokes  and  Padgett  [14]. 

For  a  detailed  survey  of  analytical  and  numerical  methods 
•Supported  by  U.S.  Army  Research  Contract  no.  DAAG29-85-G-0109 
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Abstract 

Chebyshev  method  for  solving  random  differential  equa¬ 
tion  is  presented.  The  convergence  of  the  random  coefficients 
of  the  Chebyshev  series  is  established.  Statistical  pro¬ 
perties  of  the  random  coefficients  are  discussed. 


1.  INTRODUCTION 

In  recent  years,  increasing  interest  in  the  numerical 
solution  of  random  differential  equations  has  led  to  the  pro¬ 
gressive  development  of  several  numerical  methods.  A  large 
number  of  papers  have  appeared  in  the  literature  containing 
approximate  solutions  of  random  differential  equations.  For 
Newton's  method,  successive  approximations,  perturbation 
methods,  method  of  moments,  finite  element  methods  and  other 
methods  on  the  approximate  solution  of  random  equations  we 
refer  to  Bharucha-Reid  {II.  Numerical  methods  of  random  poly¬ 
nomials  can  be  found  in  Bharucha-Reid  and  Sambandham  [2] . 

For  a  short  and  elegant  note  on  several  of  these  methods  we 
also  refer  to  Lax  [9,  10].  Some  analytical  and  numerical 
estimates  on  error  estimates  of  stochastic  differential 
equations  are  presented  in  Ladde  et  al.  [6-8].  For  other 
interesting  numerical  techniques  we  refer  to  Boyce  [3]  and 
[11] .  Numerical  treatment  of  Ito  equations  can  be  found  in 
Klauder  and  Petersen  [4],  Rumelin  [12]  and  Taley  [15].  We 
notice  that  most  of  these  numerical  methods  are  successful 

*Research  suoported  by  U.S.  Army  Research  Office,  Grant  No. 
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A  mixed  min-max  control  problem  is  considered  using 
relaxed  controls.  Mixed  min-max  problems  governed  by  dif¬ 
ferential  equations  have  been  considered  using  Dubovitskii 
and  Milyutin  Theory  and  Convex  Analysis  [6] .  Here  we  use 
relaxed  controls  and  penalization  to  deal  with  a  process 
governed  by  integral  equations.  We  use  a  technique  employed 
by  us  in  dealing  with  control  problems  governed  by  integral 
equations  [2] . 
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ABSTRACT 

Statistical  properties  of  the  numerical  solutions  of  random  difference 
equations  are  estimated.  By  an  application  of  variation  of  constant 
formula,  error  estimates  between  random  solutions  and  smooth  solutions 
(deterministic  solution)  are  discussed. 


1.  INTRODUCTION 

Numerical  solutions  of  mathematical  models  of  dynamical  systems  In  applied 
mathematics  demand  a  fair  knowledge  of  stochastic  difference  equations. 
That  is  stochastic  difference  equations  where  some  of  the  variables  can 
change  stochastically  In  time.  More  recently  Ladde  and  Sambandham  [5] 
developed  several  random  difference  Inequalities  which  are  very  powerful 
tools  to  study  stability  properties  of  stochastic  difference  equations. 

By  an  application  of  variation  of  constant  formula,  error  estimates 
between  stochastic  systems  and  the  respective  deterministic  (smooth)  sys¬ 
tems  are  discussed  In  [6].  These  estimates  provide  statistical  properties 
of  the  upper  bound  for  the  error  estimates.  For  related  results  we  refer 
to  Fal  Ma  and  Caughey  [2,3,4],  Mann  and  Wald  [7],  Deller  [1]. 

We  organize  our  article  as  follows.  In  Section  2  we  present  the 
analytical  upper  bounds  for  error  estimates.  In  Section  3  we  Include  a 
few  discussions  on  the  numerical  solutions.  In  addition,  we  have  pre¬ 
sented  a  few  tables  and  figures  to  illustrate  the  behavior  of  the  error 
between  the  mean  of  the  random  solutions  and  deterministic  solutions.  We 
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Abstract 


In  this  article  we  present  a  brief  sum¬ 
mary  of  the  numerical  solution  of  a  system 
of  random  Volterra  integral  equations.  The 
methods  we  use  are  (i)  Newton's  method  and 
(ii)  successive  approximation  method.  Based 
on  the  simulation,  we  discuss  the  mean  and 
variance  of  the  solution  of  a  system  of  random 
Volterra  integral  equations. 
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GAS  FLOWS  generate:  by  solid-propellant 

BURNING  UNDER  RANDOM  ENVIRONMENT. 


In  the  literature  it  is  customary  t  .  as¬ 
sume  that  the  local  mass  source  term, 
momentum  source  term  and  energy  source 
term  to  be  deterministic.  In  this  paper 
we  consider  some  of  these  source  terms  to 
be  random.  We  assume  that  the  flow  gener¬ 
ated  by  gasification  of  the  solid  pro¬ 
pellant  material  is  one  dimensional  and 
unsteady.  By  an  application  of  random 
choice  method,  we  discuss  the  statistical 
properties  of  all  field  variables. 
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Optimal  Processes  Governed  by  Integral  Equations 
with  Unilateral  Constraint 


Negash  G.  Medhin* 


An  optimal  process  with  unilateral  constraint  is  con¬ 
sidered  using  relaxed  controls.  Existence  of  optimal 
control  is  automatic.  Further,  problems  with  linear  and 
certain  quadratic  cost  functions  with  linear  integral  and 
state  constraints  are  such  that  the  controls  are  actually 
ordinary,  i.e.,  bounded  measurable  functions.  In  addition, 
if  the  control  set  is  a  convex  polyhedron  the  optimal  con¬ 
trols  take  their  values  in  the  vertices  of  the  polyhedron. 
Control  processes  governed  by  integral  equations  have  been 
considered  in  [2].  However,  in  [2]  existence  is  not  shown 
for  the  general  problem  and  optimality  conditions  are  not 
available  for  the  problem  when  state  constraints  are  pre¬ 
sent.  Furthermore,  the  work  in  [2]  has  some  errors  as 
pointed  out  in  [3].  Finally  we  point  out  the  basic  idea  in 
our  approach  in  this  paper- is  the  same  as  in  [6]. 
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I  A  mixed  min-max  control  problem  is  considered  using 

k 

f  relaxed  controls.  Mixed  min-max  problems  governed  by  dif¬ 

ferential  equations  have  been  considered  using  Dubovitskii 
and  Milyutin  Theory  and  Convex  Analysis  [6].  Here  we  use 
relaxed  controls  and  penalization  to  deal  with  a  process 
governed  by  integral  equations.  We  use  a  technique  employed 
[  by  us  in  dealing  with  control  problems  governed  by  integral 

equations  (2J. 
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Introduction 

We  consider  a  min-max  optimization  problem  using  relaxed 
controls.  Such  problems  have  been  considered  using 
Dubovitskii  and  Milyutin  Theory  and  convex  analysis  [5] . 

We  reformulate  the  problem  where  a  unilateral  constraint  is 
added.  We  deal  with  the  reformulated  problem  using  penali¬ 
zation  [2],  [3].  Our  procedure  produces  more  optimality 
conditions  than  [  5]  and  also  more  insight  into  (  5]  .  In  the 
special  case  where  the  controls  appear  linearly  and  the  con¬ 
trol  sets  are  convex  polyhedra  the  relaxed  controls  are 
ordinary  controls  taking  values  in  the  vertices  of  the 
polyhedra.  Finally,  we  present  examples  where  we  demonstrate 
how  the  optimal  controls  and  trajectories  could  be  calculated 
and  some  that  are  numerically  worked  out. 


KEY  WORDS:  Relaxed  controls,  unilateral  constraint, 
penalization 
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Abstract 


In  this  article  we  present  a  brief  sum¬ 
mary  of  the  numerical  solution  of  a  system 
of  random  Volterra  integral  equations.  The 
methods  we  use  are  (i)  Newton's  method  and 
(ii)  successive  approximation  method.  Based 
on  the  simulation,  we  discuss  the  mean  and 
variance  of  the  solution  of  a  system  of  random 
Volterra  integral  equations. 
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Abstract 

Chebyshev  method  for  solving  random  differential  equa¬ 
tion  is  presented.  The  convergence  of  the  random  coefficients 
of  the  Chebyshev  series  is  established.  Statistical  pro¬ 
perties  of  the  random  coefficients  are  discussed. 


1.  INTRODUCTION 

In  recent  years,  increasing  interest  in  the  numerical 
solution  of  random  differential  equations  has  led  to  the  pro¬ 
gressive  development  of  several  numerical  methods.  A  large 
number  of  papers  have  appeared  in  the  literature  containing 
approximate  solutions  of  random  differential  equations.  For 
Newton's  method,  successive  approximations,  perturbation 
methods,  method  of  moments,  finite  element  methods  and  other 
methods  on  the  approximate  solution  of  random  equations  we 
refer  to  Bharucha-Reid  (1).  Numerical  methods  of  random  poly¬ 
nomials  can  be  found  in  Bharucha-Reid  and  Sambandham  [2]. 

For  a  short  and  elegant  note  on  several  of  these  methods  we 
also  refer  to  Lax  (9,  10].  Some  analytical  and  numerical 
estimates  on  error  estimates  of  stochastic  differential 
equations  are  presented  in  Ladde  et  al.  [6-8] .  For  other 
interesting  numerical  techniques  we  refer  to  Boyce  [3]  and 
[11] .  Numerical  treatment  of  Ito  equations  can  be  found  in 
Klauder  and  Petersen  [4],  Rumelin  [12]  and  Taley  [15].  We 
notice  that  most  of  these  numerioal  methods  are  successful 
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ABSTRACT 

The  solution  of  several  elasticity  problems,  and  parti¬ 
cularly  crack  problems,  can  be  reduced  to  the  solution  of 
one-dimensional  singular  integral  equations  with  a  Cauchy- 
type  kernel.  In  this  paper,  we  present  a  method  for  the 
numerical  solution  of  a  random  singular  integral  equation 
of  the  Cauchy-type.  To  illustrate  this  method,  it  is  applied 
to  the  random  singular  integral  equation  that  arises  in  the 
problem  of  periodic  array  of  straight  cracks  inside  an 
infinite  isotropic  elastic  medium  and  subjected  to  random 
pressure  distribution  along  the  crack  edges.  The  statistical 
properties  of  the  random  solution  are  evaluated  numerically, 
and  used  to  determine  the  stress  intensity  factor  at  the 
crack  tips.  Results  of  this  study  highlight  the  advantage  of 
using  this  method  to  solve  a  random  Cauchy-type  singular 
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integral  equation.  The  method  is  applied  to  the  crack  problem 
subjected  to  different  forcing  functions. 


1.  INTRODUCTION 

In  recent  years,  an  unprecedented  widespread  interest 
in  the  numerical  solution  of  singular  integral  equations  of 
the  Cauchy-type  has  resulted  in  the  development  of  several 
numerical  methods  [1-5]  besides  stimulating  a  mathematical 
interest  in  their  numerical  analysis.  Since  the  pioneering 
work  by  Erdogan  [1]  in  1969  and  subsequently  by  Erdogan  and 
co-workers  [3],  a  large  number  of  papers  have  appeared  in  the 
literature  concerning  singular  integral  equations  that  arise 
in  several  fields  of  mathematical  physics  and  engineering. 

The  solution  methods  developed  have  aimed  at  reducing  the 
singular  integral  equation  to  a  system  of  linear  algebraic 
equations  which  can  be  easily  solved  to  give  an  approximate 
expression  for  the  solution  of  a  singular  integral  equation  or 
solutions  in  the  case  of  a  system  of  singular  integral  equa¬ 
tions  . 

Several  of  the  numerical  methods  developed  over  the  years 
were  found  to  have  disadvantages  due  to  their  inherent  com¬ 
plexity,  low  degree  of  accuracy  and  limitations  as  regards  the 
classes  of  singular  integral  equations  to  which  they  were 
applicable.  It  is  known  that  numerical  solutions  of  Fredholm 
integral  equations  by  reduction  to  a  system  of  linear  equations 
could  be  extended  to  the  case  of  Cauchy  singular  integral 


\  ‘  . 


3 


equations  (CSIEs),  that  is  to  the  case  when  the  kernels  of 
the  Fredholm  integral  equations  have  Cauchy-type  singularities. 
Since  Cauchy-type  singular  integral  equations  on  closed  con¬ 
tours  have  been  treated  elsewhere  [7,  8],  in  this  paper  we 
confine  our  attention  to  equations  on  intervals  (a,b)  of  the 
real  line,  which  may  be  either  finite  or  infinite.  In  parti¬ 
cular,  we  consider  linear  equations  since  the  theory  for  these 
is  fairly  complete  and  well  developed.  The  most  frequently 
used  technique  for  the  numerical  solution  of  linear  singular 
integral  equations  consists  in  approximating  the  integral 
terms  by  using  an  approximate  numerical  integration  rule 
and  applying  the  singular  integral  equation  at  appropriately 
selected  collocation  points  (2-4,  9].  One  determination  of 
the  unknown  function  at  the  abscissas  used,  the  function  is 
approximated  along  the  whole  integration  interval  [10  ]. 

A  wide  body  of  literature  pertaining  to  the  numerical 
solution  of  singular  integral  equations  of  the  deterministic 
type  is  succintly  summarized  by  Ergodan,  et  al.  [3],  Golberg 
[11],  and  Theocaris  [6].  In  referring  to  singular  integral 
equations,  Erdogan,  Golberg  and  Theocaris  in  their  articles 
imply  Cauchy-type  singular  integrals  and  Cauchy-type  singular 
integral  equations.  We  refer  to  Bharucha-Reid  and  Bharucha- 
Reid  and  Christensen  [12,  13]  for  surveys  of  the  different 
analytical  methods  for  the  solution  of  random  integral  equa¬ 
tions.  Most  recently,  Christensen  and  Bharucha-Reid  [14] 
developed  numerical  solution  procedures  for  Fredholm  equations 
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with  random  kernels,  random  right  hand  side  or  both  and  for 
Fredholm  equations  with  random  degenerate  kernels  [15].  The 
Cauchy-type  singular  integral  equations  which  arise  fre¬ 
quently  in  the  fields  of  aerodynamics,  hydrodynamics,  solid 
mechanics,  elasticity  and  fracture  mechanics,  are  closely 
related  to  Fredholm-type  equations.  There  exists  a  pragmatic 
need  to  develop  numerical  techniques  in  order  to  be  able  to 
study  the  characteristics  of  the  random  singular  integral 
equations.  The  rationale  for  this  study  was  to  extend  the 
concept  for  the  solution  of  random  Fredholm  integral  equations 
[14-16],  to  the  solution  of  random  singular  integral  equations 
of  the  Cauchy-type  that  arise  in  crack  problems  in  the 
classical  theory  of  elasticity.  In  particular,  we  consider 
the  problem  of  an  array  of  periodic  cracks  inside  an  infinite 
isotropic  elastic  medium  and  subjected  to  random  pressure 
distribution  along  the  crack  edges,  i.e.,  the  case  of  a 
random  forcing  function.  The  statistical  properties  of  the 
solution  are  evaluated  and  the  results  obtained  for  the  random 
equation  are  compared  with  those  of  the  deterministic  equa¬ 
tion.  Further  we  also  discuss  the  distribution  of  the  stress 
intensity  factor  for  different  sample  sizes. 

2.  SINGULAR  INTEGRAL  EQUATIONS 

Consider  the  singular  integral  equation  of  the  form 

(2.1)  i  ;  dt  ♦  k(x,t)  }(t)dt  .  f  (x>  ,  -1  <  X  <  1, 

’-it-*  -1 
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where  k(x,t)  is  a  regular  kernel.  (It)  is  the  unknown  func¬ 
tion  proportional  to  the  crack  tip  opening  displacement  and 
f (x)  is  the  known  function  representing  pressure  distribu¬ 
tion  along  the  crack. 

Let  g(t)  =  4)(t)/z(t),  where  z(t)  =  (1-t^)  and 


(2.2) 


/  <i>(t)dt  =  0. 
-1 


We  consider  (2.1)  together  with  (2.2)  and  replace  g(t)  by  a 
new  unknown  function  G(x)  defined  by  the  first  integral  of 
(2.1) ,  that  is , 


(X)  =  4  /  (l-t^)'^-^^  dx. 


-1  <  X  <  1 . 


Since  g(t)  satisfies  (2.2),  we  can  express  g(t)  in  terms  of 


G(x)  by 


(2.4)  g(t)  =  T  ■'  (1-y^)^^^ 


t  -y 


dy,  -1  <  t  <  1. 


By  an  application  of  (2.3)  and  (2.4),  we  can  write  (2.1)  in 
the  form 

^  2  1/2 

(2.5)  G(x)  +  /  (1-y  )  .  K(y,x)G(y)dy  =  f (x)  -1  <  x  <  1, 


where 


(2.6) 


K(y,x)  =  -  /  (1-t^)"^/^  dt 

-1  ^  ^ 


We  remark  that  (2.5)  is  another  form  of  (2.1).  The  equation 
(2.5)  is  Fredholm  equation  of  second  type  and  its  numerical 
solution  can  be  found  by  quadrature  method. 
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The  probabilistic  analogue  of  (2.5)  is  the  random  equa¬ 
tion.  We  call  an  equation  random  if  certain  components  such 
as  coefficients,  kernels,  nonhomogeneous  terms  or  forcing 
functions  and  initial  and  or  boundary  data  are  random  func¬ 
tions.  The  random  version  of  (2.5)  is  of  the  form 


(2.7)  G(x,lo)  +  /  (l-y^)  ^/^K(y,x,uj)G(x,u))dy  =  f(x,to) 

-1 

-1  <  x  <  1 

where  K(y,x,w)  is  the  random  kernel  and  f{x,uj)  is  the  random 
forcing  term.  The  parameter  oj  is  an  element  of  a  given 
probability  space  (Q,A,P).  Here  G(x,u))  is  the  random  func¬ 
tion  referred  to  as  the  solution  or  output.  For  numerical 
solution  of  (2.7)  we  refer  to  (14,  15,  16]. 

In  this  article  we  consider  a  periodic  array  of  cracks 
along  a  straight  line  in  an  infinite  isotropic  elastic  medium 
under  plane  strain  or  generalized  plane  stress.  The  length 
of  each  crack  is  assumed  to  be  '2a'  and  the  period  of  the 
array  is  'b'.  For  t^'e  case  of  a  constant  compressive  loading 
distribution  along  the  crack  edges,  loakimidis  [17]  found  the 
singular  integral  equation  to  be 


(2.8) 


a 

5 


(1-t 


2-1/2 


cot 


-a (t-x) 


)  g  (t) dt  =  f (x)  -1 


-1 


which  is  of  the  form  (2.1)  with 


(2.9) 


,  , ^  ,  a  ,"a (t-x) ^  1 

k(t,x)  =  5  cot( - g - )  - 


where  cracks  are  assumed  to  be  loaded  by  constant  compressive 
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loading  distribution  f(x)  along  both  crack  edges.  If  the 
loading  distribution  f (x)  is  random,  then  (2.8)  is  a  random 
singular  integral  equation.  If  we  denote  random  loading 
distribution  by  f(x,uj)  then  the  random  version  of  (2.8)  is 

(2.10)  g  /  ^  f(x,aj)  -1  <  x  <  1 

Equation  (2.10)  together  with  (2.2)  can  be  solved  numerically 
by  methods  suggested  in  [12-17]. 


3.  NUMERICAL  TREATMENT 

In  this  section  we  discuss  the  numerical  solution  of 
(2.10) .  We  assume  the  ratio  a/b  =  .4  and  the  successive 
cracks  lie  close  to  each  other.  We  denote  the  value  of  the 
stress  intensity  factor  by  K(±l,uj)  at  the  crack  tips  t  =  ±1. 
It  is  given  by 


(3.1)  K(2:1,u))  =  ±g(±l,u)) 

1  ^ 

~  N  \  ±x^)G(x^,u.)  . 

k=l 

where  x,  are  suitable  nodes 
k 

Xj^  =  cos(k-0.5)  ^ 

which  we  use  to  solve  (2.7).  The  unknown  random  function 
g(x,uj)  in  (2.10)  is  evaluated  at  different  pressure  distri¬ 
butions  along  the  crack  edge?,  namely. 


Example  3.1. 


f  (X,ij)  e  N  (1 ,0^) 
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f(x,^)  £  M(e^,a^) 
f(x,^)  £  N(e^,a^e^^) 


A  sample  of  5000  standard  normal  random  numbers  are 
generated  from  IMSL  routines  (GGNML) .  The  resulting  random 
numbers  are  used  to  generate  the  random  input  functions, 
namely,  the  random  pressure  distributions  in  the  manner 
described  above.  For  each  random  input  function,  equation 
(2.7)  and  hence  equation  (2.10)  is  solved  numerically  to 
obtain  the  random  output  function,  and  the  stress  intensity 
factor.  At  the  crack  tips,  K(il,(ju)  are  computed  based  on  the 
relationship  in  equation  (3.1). 


4.  RESULTS  AND  DISCUSSION 

Accurate  estimates  of  the  stress  intensity  factor  at  the 
crack  tips  under  more  realistic  pressure  distributions  is 
essential  for  the  reliable  prediction  of  fatigue  crack  growth 
rates  in  structural  members.  Several  analytical  techniques 
are  presently  available  to  obtain  solutions  to  random 
Cauchy-type  singular  integral  equations  that  arise  in  crack 
problems.  Most  of  these  numerical  techniques  for  obtaining 
the  stress  intensity  factors  for  crack  and  elasticity  prob¬ 
lems  have  been  concerned  with  pressure  distributions  that  are 
purely  deterministic  in  nature. 

In  dealing  with  random  equations,  the  main  objectives 
to  be  concerned  with  are:  (1)  determination  of  the  statisti¬ 
cal  properties  of  the  random  solution  such  as  its  expectation. 
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variance  and  higher  moments;  (2)  to  establish  the  relation¬ 
ship  between  the  expected  solution  of  the  random  equation 
and  the  solution  of  the  deterministic  equation  or  the  mean 
solution  [18];  (3)  determining  the  distribution  of  the  solu¬ 

tions,  and  (4)  discussing  the  limiting  properties  of  the 
solutions.  Of  particular  interest  is  the  relationship 
between  the  statistical  properties  of  the  solution  and  the 
statistical  properties  of  the  random  input  function  intro¬ 
duced  into  the  equation.  Table  I-V  summarize  the  statisti¬ 
cal  properties  of  K(±l,(jo). 

The  distribution  of  stress  intensity  factor  K(il,Lu) 
for  the  functions  f(x,u)  are  shown  in  Figures  1-30.  These 
figures  show  that  as  the  value  of  N  increases,  the  distribu¬ 
tion  of  K(il,Gi))  is  approximately  normal  with  E[K(±l,(jj)]  being 
the  solution  of  the  mean  equation.  It  is  observed  that  as 
the  value  of  N  increases,  the  distribution  of  the  stress 
intensity  factor  K(±l,oj)  is  normal  as  long  as  the  forcing 
function  f(x,ji))  is  normal. 

5.  CONCLUSIONS 

1.  The  foregoing  analysis  demonstrates  a  fairly  simple  yet 
useful  method  for  the  solution  of  random  Cauchy-type 
singular  integral  equation  that  arises  in  a  crack  prob¬ 
lem  in  the  classical  theory  of  elasticity.  The  solu¬ 
tion  of  the  random  equation  helps  in  the  accurate 
evaluation  of  the  stress  intensity  factor  at  the  crack 


tips  under  more  realistic  input  or  forcing  functions. 

The  distributions  of  the  stress  intensity  factors 

are  approximately  normal  when  the  forcing  func¬ 
tions  are  normal. 

The  results  of  this  study  highlight  the  influence  of  the 
nature  of  the  random  input  function  on  the  stress 
intensity  factor. 


TABLE  I 


f  (  X  ,  - ) 

e  N(l,0.02^) 

Mean  of 

Variance  of 

Sample  Size 

N 

K(l,aj) 

K  (l,uj) 

500 

2 

1.634664 

0.001065 

3 

1.575876 

0.000896 

5 

1. 563470 

0.000866 

7 

1.566512 

0.000969 

9 

1.566098 

0.000910 

1000 

2 

1.635569 

0.001079 

3 

1.575410 

0.000960 

5 

1.  563799 

0.000931 

7 

1.567005 

0.001003 

9 

1.565704 

0.000940 

2000 

2 

1.636609 

0.001098 

3 

1.575890 

0.000962 

5 

1.564721 

0.000965 

7 

1.566515 

0.000993 

9 

1. 564965 

0.000987 

3000 

2 

1.637003 

0.001088 

3 

1.575030 

0.000962 

5 

1.564731 

0.000977 

7 

1.566254 

0.000990 

9 

1.564981 

0.000974 

4000 

2 

1.636695 

0.001094 

3 

1.575152 

0.000984 

5 

1. 565015 

0.000977 

7 

1.566710 

0.000985 

9 

1.564917 

0.000983 

5000 

2 

1.636726 

0.001089 

3 

1.575127 

0.000969 

5 

1.564924 

0.000989 

7 

1. 566634 

0.000969 

9 

1.564916 

0.000983 

■  »  “■--  *.*  "j*  “■/•  '■• 
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TABLE  II 


f{x,^)  £  N(e  ,0.02") 


Mean  of 

Variance  of 

Sample  Size 

N 

K'l,.j) 

K{l,u)) 

500 

2 

2.567233 

0.000985 

3 

2.548983 

0.000967 

5 

2.538348 

0.000880 

7 

2.539721 

0.000997 

9 

2.539926 

0.001007 

1000 

2 

2.566281 

0.001026 

3 

2.548836 

0.001010 

5 

2.537966 

0.000897 

7 

2.540643 

0.000933 

9 

2.539931 

0.000971 

2000 

2 

2.566361 

0.001063 

3 

2.548886 

0.001033 

5 

2.538958 

0.000955 

7 

2.540021 

0.000937 

0 

2.539534 

0.000992 

3000 

2 

2.566199 

0.001082 

3 

2.549063 

0.001023 

5 

2.539601 

0.000987 

7 

2.539973 

0.000944 

9 

2.539349 

0.000992 

4000 

2 

2.565732 

0.001077 

3 

2.548748 

0.001026 

5 

2.539587 

0.000976 

7 

2.539796 

0.000949 

9 

2.539454 

0.001007 

5000 

2 

2.565745 

0.001069 

3 

2.548838 

0.001015 

5 

2.539710 

0.000975 

7 

2.539714 

0.000955 

9 

2.539338 

0.001008 

TABLE  III 


f  (X,uj) 

X  2  ^  V 

e  N(e  ,0.02  e 

Mean  of 

Variance  of 

Sample  Size 

N 

K(l,ai) 

K(l,aj) 

500 

2 

2.565219 

0.002773 

3 

2.544691 

0.002617 

5 

2.539554 

0.002960 

7 

2.541145 

0.002661 

7 

2.542362 

0.002784 

1000 

2 

2.565728 

0.002727 

3 

2.546564 

0.002680 

5 

2.541254 

0.002728 

7 

2.541461 

0.002704 

9 

2.542656 

0.002535 

2000 

2 

2.565767 

0.002784 

3 

2.548370 

0.002618 

5 

2.539482 

0.002647 

7 

2.540158 

0.002643 

9 

2.540600 

0.002528 

3000 

2 

2.565145 

0.002694 

3 

2.547493 

0.002589 

5 

2.539537 

0.002630 

7 

2.539658 

0.002603 

9 

2.541302 

0.002609 

4000 

2 

2.565149 

0.002677 

3 

2.547726 

0.002611 

5 

2.539973 

0.002641 

7 

2.538959 

0.002602 

9 

2.540734 

0.002614 

5000 

2 

2.564995 

0.002644 

3 

2.547745 

0.002592 

5 

2.539966 

0.002628 

7 

2.539020 

0.002585 

9 

2.540420 

0.002630 

TABLE  IV 


f  (X,uj) 

£  N(e~^,0. 

2 

02^) 

Mean  of 

Variance  of 

Sample  Size 

N 

K(-l,uj) 

K(~l,to) 

500 

2 

1.562386 

0.000996 

3 

1.478523 

0.001107 

5 

1.465161 

0.000979 

7 

1.466776 

0.001004 

9 

1.469071 

0.000901 

1000 

2 

1.563262 

0.001020 

3 

1.479762 

0.001065 

5 

1.465880 

0.001015 

7 

1.466557 

0.000990 

9 

1.468161 

0.000939 

2000 

2 

1.563492 

0.001015 

3 

1.479095 

0.000998 

5 

1.466543 

0.001003 

7 

1.466125 

0.000990 

9 

1.468237 

0.000972 

3000 

2 

1.562768 

0.001029 

3 

1.478868 

0.000993 

5 

1.467130 

0.000992 

7 

1.466727 

0.000992 

9 

1.468084 

0.000990 

4000 

2 

1.562527 

0.001035 

3 

1.478661 

0.000985 

5 

1.467698 

0.000993 

7 

1.466606 

0.000992 

9 

1.467634 

0.000986 

5000 

2 

1.562371 

0.001046 

3 

1.478330 

0.000980 

5 

1.467851 

0.000987 

7 

1.466590 

0.000985 

9 

1.467648 

0.000978 

TTTV 
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TABLE  V 


f  (x, 

.^)  - 

'  n 

N(e‘^,0.02" 

r^T3r 

e  ) 

Mean  of 

Variance  of 

Sample  Size 

N 

K(-l,u)) 

K(-l,a)) 

500 

2 

1.563375 

0,000957 

3 

1.480107 

0.000776 

5 

1.468341 

0.000749 

7 

1.467878 

0.000831 

9 

1.467650 

0.000818 

1000 

2 

1.562688 

0.000943 

3 

1.479829 

0.000833 

5 

1.467521 

0.000833 

7 

1.466508 

0.000835 

9 

1.468472 

0.000826 

2000 

2 

1.561791 

0.000963 

3 

1.479141 

0.000852 

5 

1.467389 

0.000858 

7 

1.467030 

0.000867 

9 

1.467370 

0.000843 

3000 

2 

1.561813 

0.000973 

3 

1,478927 

0.000861 

5 

1.467764 

0.000881 

7 

1.466812 

0.000870 

9 

1.467162 

0.000861 

4000 

2 

1.562072 

0.000976 

3 

1. 478724 

0.000376 

5 

1.467863 

0.000880 

7 

1.467345 

0.000872 

9 

1.467029 

0.000851 

5000 

' 

1.561947 

0.000975 

1.478948 

0.000869 

1.467711 

0.000880 

1.467042 

0.000875 

1.466858 

0.000842 

;-;i 


^■:v 

r..s 


SIZE  500 


N-7  SAMPLE  SIZE  500 


f(x,uj)  £  N{1,0.02‘‘) 


1.65 

|1,  id) 
g,  3. 


PLE  SIZE  500 


f(x,u.)  £  N(l,0.02‘) 


1.65 
1,  id) 


flx,^)  t  N(l,0.02-‘) 


N-7  SAMPLE  SIZE  5000 


3 

O 

i.i 

S  0  1 

Ll 


f(x,--)  €  N(l,0.02^) 


£  0.1 
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ABSTRACT 


In  this  paper,  we  present  a  aethod  for  the  nuaerlcal  solution  of  randoa 
Lowe's  Integral  equation.  The  technique  Involves  a  Chebyshev  series  approxlea- 
tlon,  the  coefficients  of  which  are  the  solutions  of  a  system  of  random  linear 
equations.  The  statistical  properties  of  the  random  solution  are  evaluated 
nuaerlcally  and  used  to  determine  the  distribution  of  the  random  coefficients. 
We  express  these  distributions  in  histograms. 


K,"  H.' 


'  "w  ’S.^  “-*»«  "i,"^ 


“.  ^  N.’  V  -..^ 


Abstracts:  Contribated  Papers  and  Poster  Presentations 


algebras  L  a  line  of  t.  The  combvnatorial 

properties  of  the  (0.1)  -  matrices  which 

span  these  algebras  are  transferred  to  a 
polynomial  ring  R  =  kCx^ , . . . ,x^)/I .  I  a  certain 

homogeneous  ideal.  R  is  noetherian  of  Krull 
dimension  n* 1  in  which  certain  statements  about 
the  prime  ideals  of  R  are  equivalent  to 
statements  about  the  co 11 ineations  of  The 
Krull  intersection  theorem  provides  an  interesting 
topology  for  R  when  i  admits  no  col lineation. 

John  L.  Hayden 

Bowling  Green  State  University 
Department  of  Mathematics 
Bowling  Green,  OH  43403 


«38  (C/P  Session  1;  Wed.  12:OONOON) 
Family  of  Algorithms  for  the  Graph  Bisect 
Problem 


The  Graph  Bisection  Problem,  GBP,  is  the  problem 
of  partitioning  the  vertex  sec  of  a  given  graph 
into  two  subsets  of  Che  same  sice  (within  one 
element)  so  Chat  the  number  of  edges  cut  is 
minimi  zed . 

We  present  a  family  of  algorithms  for  GBP  that 
implements  Che  idea  of  a  restricted  branch  and 
bound  method. 

I:  is  shown  that  for  every  constant  \>0,  there  is 
an  algorithm  with  running  time  of  0(c")  (c<2), 
which  performs  as  an  exact  algorithm  on  the  sec  of 
graphs  with  n  vertices  and  Sin  edges. 

We  present  the  results  of  computer  simulations 
with  the  algorithms  and  discuss  several  conjec¬ 
tures  which  are  based  on  these  results. 


Mark  R.  Goldberg 
Department  of  Computer  Science 
Rensselaer  Polytechnic  Institute 
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<39  (C/?  Session  8;  Thu.  12:OONOON) 

APPROX  IM.ATIO.S  TO  AuTOREGRESIVE  .MODEL  WITH 
STOCHASTIC  COEFFICIENTS 

An  autoregressive  (.4R)  model  with  stochas¬ 
tic  cofficients  is  generally  modelled  by 
deterministic  AR  node!  where  stochastic 
coefficients  are  replaced  by  the  mean  of 
the  stochastic  coefficients.  It  has  been 
proved  that  this  approximation  is  restric¬ 
tive  and  not  applicable  in  general.  In 
this  article  we  suggest  suitable  improved 
appr  oxijiat  i  ons  which  are  less  restrictive 
and  more  'closer'  to  the  exact  value.  Some 
applications  of  pcwer  spectral  estimation 
algorithm  with  stochastic  coefficients  are 


also  discussed.  We  illustrate  our  analy¬ 
tical  techniques  through  computer  algori¬ 
thms  and  numerical  results.  These  results 
will  highlight  the  error  involved  if  one 
uses  deterministic  .AR  models  insted  of 
stochastic  .AR  models. 

M.S.AMBANDHAM 

Department  of  Mathematics 

Morehouse  College/  Atlanta  University 

Atlanta  GA  30314 


K.B.  BOTA 

Department  of  Physics 
Atlanta  University 
Atlanta  GA  30314 


#90  (C/P  Session  13;  Fri .  5:0QPM) 
Discriminative  Optimization 

The  excentions  to  Che  simplex  algorithm 
required  by  Goal  Programming  (GP)  in¬ 
spire  a  further  generalization  on  a 
broader  class  of  mixed  binary  problems 
with  a  structure  among  goal  functions. 

The  Discriminative  Optimization  (DO) 
is  defined  in  it's  generality  and  the 
score  of  optimization  methods  like;  Set¬ 
up  cost,  Goal  Tree  Programming,  Aggrega¬ 
tive  Goal  Programming,  Linear  Goal  Prog¬ 
ramming,  etc.  are  shown  Co  be  particular 
cases  of  Che  same  scheme.  Correspondig- 
ly,  simplex-based  solution  algorithms 
for  Chose  problems,  and  even  for  Che 
Linear  Programming,  can  be  obtained  as 
implementations  of  the  general  solution 
process . 

Dr.  Pawel  Radzikowski 
Seton  Hall  University 
School  of  Business 
South  Orange,  NJ  07079 


*91  (C/P  Session  15  ;  Fri.  5:0QPM) 

A  Model  for  Cannibalism 


Many  species  of  insects  and  fish  practice 
cannibalism,  presumably  due  to  selective 
pressure  toward  this  behavior.  This  paper 
pp^sents  a  system  of  non-linear  ODE’s  which 
reflects  population  dynamics  in  an  insect 
species  whose  members  pass  through  a  sequence 
of  stages,  and  in  which  later  stages 
cannibalize  earlier  stages.  Solution 
trajectories  must  converge  when  the  number 
of  stages  is  four  or  loss.  conjec-ure 
that  convergence  always  -obtains  but,  wnen 
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ABSTRACT 


APPLIED  MATHEMATICS 

SEKO,  JOSEPH  T.  M.S .  ATLANTA  UNIVERSITY,  1986 

THE  USE  OF  LAGRANGE  MULTIPLIERS  AND  KUHN-TUCKER* S 
THEORY  IN  OPTIMIZATION  PROBLEMS 

Advisor:  Professor  Negash  Medhin 

Thesis  dated  November,  1986 

Lagrange  multipliers,  penalty  methods,  and  Kuhn- 
Tucker's  theory  are  some  important  mathematical  tools  used 
in  optimization  problems.  These  tools  are  discussed  so 
that  one  can  appreciate  the  current  areas  of  optimization 
research.  Moreover,  since  extensive  research  work  exists 
for  linear  optimization  problems,  only  nonlinear 
applications  are  discussed. 
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EPHRAIM,  DANIEL  B.S,,  MOREHOUSE  COLLEGE,  1985 

NUMERICAL  SOLUTION  OF  LINEAR  INTEGRAL  EQUATIONS 
WITH  RANDOM  FORCING  TERMS 

Advisor:  M.  Sambandhatn 

Thesis  dated  July,  1987 

In  Chapter  one  of  this  report  we  define  Fredholm  integral 
equations  of  the  second  kind,  Volterra  integral  equations  of 
the  second  kind  and  differentiate  between  the  two  of  them 
and  explain  why  integral  equations  are  important.  In  Chapter 
two  we  discuss  numerical  procedures  to  integral  equations. 

The  equations  we  used  in  this  report  are  of  two  types:  (1) 
Fredholm  equations  and  (2)  Volterra  equations.  The  methods 
we  used  for  Fredholm  equations  are:  (i)  Simpson's  rule,  (ii) 
Trapezoidal  rule,  (iii)  Weddle's  rule,  (iv)  the  Collocation 
method,  and  (iv)  the  Galerkin  method. 

For  Volterra  equations  we  used  the  successive  approxima¬ 
tion  method  with  (i)  Simpson's  rule,  (ii)  Trapezoidal  rule 
and  (iii)  Weddle's  rule  to  evaluate  the  integrals. 

In  both  Fredholm  and  Volterra  integral  equations  we  have 
the  forcing  term  to  be  random.  Our  simulation  results  are 
presented  in  tables  and  graphs. 


